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ABSTRACT 

Kilohertz quasi-periodic X-ray brightness oscillations (kilohertz QPOs) have now been de- 
tected in more than twenty accreting neutron stars in low-mass binary systems. Two kilohertz 
QPOs are usually detected in each star. Burst oscillations and two kilohertz QPOs have re- 
cently been detected in the 401 Hz accretion-powered X-ray pulsar SAX J1808.4— 3658. In this 
star the frequency of the burst oscillation is approximately equal to the star's spin frequency 
i^spin whereas the frequency separation of the two kilohertz QPOs is approximately i'spin/2. If 
as expected the frequency of the burst oscillations in other stars is also approximately i^spin; 

the 

frequency separation is w i^spin in some stars but « i/spin/2 in others. A frequency separation 
approximately equal to I'spin/S is unexplained in all existing models of the kilohertz QPOs. Here 
we propose a modified version of the sonic-point beat-frequency model that can explain within a 
single framework why the frequency separation is close to Vgpin in some stars but close to ^'spin/2 
in others. As in the original sonic-point model, the frequency vqpo2 of the upper kilohertz QPO 
is close to the orbital frequency Voih the radius r^p of the sonic point in the disk flow. We 
show that magnetic and radiation fields rotating with the star will preferentially excite vertical 
motions in the disk at the "spin-resonance" radius r^r where forb — t'spin is equal to the vertical 
epicyclic frequency. If the flow at rsr is relatively smooth, the vertical motions excited at rsr 
modulate the X-ray flux at fQpoi ~ '^qpo2 — I'spin- If instead the gas at r^r is highly clumped, 
the vertical motions excited at rsr modulate the X-ray flux at i^qpoi ~ '^qpo2 — '^spin/2. This 
sonic-point and spin-resonance model can also explain quantitatively the decrease of the kilohertz 
QPO frequency separation with increasing accretion rate that is observed in many sources. 

Subject headings: accretion — relativity — stars: neutron — X-rays: binaries — X-rays: bursts 
— X-rays: stars 
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1. Introduction 

Recent discoveries made with RXTE have clarified the relationship between the pairs of kilohertz quasi- 
periodic oscillations (kilohertz QPOs) observed in the accretion-powered X-ray emission of neutron stars in 
low-mass binary systems (see van der Klis 2000; Lamb 2002), the high-frequency oscillations seen during 
thermonuclear X-ray bursts (see Strohmayer 2001; Strohmayer & Bildsten 2003), and the spin frequency of 
the star. Key developments include the discovery of burst oscillations and kilohertz QPOs in the 401 Hz 
accretion-powered X-ray pulsar SAX J1808.4— 3658 (Chakrabarty et al. 2003; Wijnands et al. 2003) and 
burst oscillations in the 314 Hz accretion-powered X-ray pulsar XTE J1814— 338 (Markwardt & Swank 
2003; Markwardt, Strohmayer, & Swank 2003b). 

SAX J1808.4— 3658 was discovered and identified as an X-ray burst source during its September 1996 
outburst (in 't Zand et al. 1998). During its April 1998 outburst, a 401 Hz periodic oscillation ("pulsation") 
was detected in its persistent emission (Wijnands & van der Klis 1998; Chakrabarty & Morgan 1998), but 
neither kilohertz QPOs nor burst oscillations were detected then. Subsequent analysis of an X-ray burst 
recorded by BeppoSAX during the September 1996 outburst yielded a marginal detection of an oscillation 
at 400 ± 2 Hz (in 't Zand et al. 2001). SAX J1808.4-3658 was extensively observed throughout its October 
2002 outburst and periodic oscillations, burst oscillations, and kilohertz QPOs were all detected with high 
confidence (Chakrabarty et al. 2003; Wijnands et al. 2003), making SAX J1808.4— 3658 a Rosetta stone for 
understanding these phenomena. 

The results of Chakrabarty et al. (2003) establish conclusively that the spin frequency of the neutron 
star in the SAX J1808.4— 3658 system is 401 Hz: no oscillation at half this frequency was detected, and the 

upper limit on the rms amplitude of any such oscillation was determined to be 0.014%, 300 times smaller than 
the amplitude of the 401 Hz oscillation. This excludes the possibility that the spin frequency is 200.5 Hz. 

Wijnands et al. (2003) discovered a pair of kilohertz QPOs in SAX J1808.4— 3658. Their frequency 
separation Az/qpq is within 3% of i'spin/2, demonstrating that Az/qpq is related to the spin of the star and 
that the system somehow generates a frequency difference approximately equal to half the spin frequency, a 
result unexplained by existing models. 

Chakrabarty et al. (2003) observed four X-ray bursts during the October 2002 outburst of 
SAX J1808.4— 3658. All showed a ^^401 Hz oscillation, demonstrating that in this star, i^burst ~ '^spin- 
Markwardt et al. (2003b) have analyzed three X-ray bursts from XTE J1814— 338. All showed nearly coher- 
ent ~ 314 Hz oscillations, demonstrating that Uhurst ~ Jospin in this star, also. These results make compelling 
the previously strong evidence (Strohmayer et al. 1996, 1997, 1998a, 1998b; Strohmayer & Markwardt 1999; 
Muno et al. 2000; Muno et al. 2002; Strohmayer & Markwardt 2002) that burst oscillations are generated by 
the spin of the star and suggest that i^burst is close to Vspm in all the neutron stars in which burst oscillations 
have been observed. If so, then Az/qpq is f« fspin in the stars with i/gpin < 360 Hz but is ~ i'spiii/2 in the 
stars with ;/gpin > 400 Hz. A successful model must therefore explain why some stars have Az/qpo ~ Jospin, 
whereas others have Az^qpo ~ t'spin/2. 

Here we propose a modification of the original sonic-point beat-frequency model (Miller, Lamb, & 
Psaltis, hereafter MLP98; Lamb & Miller 2001) that can explain within a single framework why Afgpo is 
close to i/gpin for some stars but close to fspin/S for others. In § 2 we discuss the general implications of the 
new observations for kilohertz QPO mechanisms and in § 3 we propose specific mechanisms for generating 
the upper and lower kilohertz QPOs. We discuss our conclusions in § 4. 
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2. General Inferences from the Observations 

Taken together, the new and previous observations provide strong hints about the mechanisms that 
generate the kilohertz QPO pairs: 

1. It appears highly likely that the frequency of one of the two kilohertz QPOs reflects the orbital 

frequency of gas in the inner disk. The kilohertz QPOs have frequencies similar to those of orbital motion 
near neutron stars. The frequencies of the kilohertz QPOs vary by hundreds of Hertz on time scales as short 
as minutes (see, e.g., Mendez et al. 1999; van der Klis 2000). This is possible if their frequencies are related to 
orbital motion at a radius that varies (Lamb 2002). These frequency variations exclude mechanisms like the 
one proposed by Abramowicz & Kluzniak (2001) to explain the high-frequency QPOs observed in black hole 
candidates and by Kluzniak et al. (2003) to explain the kilohertz QPOs observed in SAX J1808.4— 3658, 
because this type of mechanism requires a low-order resonance between the geodesic frequencies of test 
particles orbiting at a fixed radius and the resonance disappears when the two frequencies change. 

2. The star's spin is somehow involved in producing the frequency separation of the two kilohertz 
QPOs in a pair. This involvement is clear in SAX J1808.4— 3658, where Ai^qpq ~ i'spin/2. It is strongly 
indicated in the other kilohertz QPO sources because in all cases where both Ai^qpo and i4>urst have been 
measured, Ai/qpq is approximately commensurable with the stellar spin frequency inferred from t^urst (see 
van der Klis 2000; Lamb 2002). In four cases, Ai/qpo is consistent with i^spin or with half Jospin to better 
than 5% (4U 1608-52: ly^pin = 319 Hz, A//qpo = 301.3 ± 7.9 Hz [Mendez et al. 1997; Mendez 2000, private 
communication]; 4U 1702—429: i^spin = 329 Hz, Ai^qpo = 333 ± 5 Hz [Markwardt, Strohmayer, & Swank 
1999]; KS 1731-260: u^pin = 524 Hz, Ai/qpq = 260 ± 10 Hz [Smith, Morgan, & Bradt 1997; Wijnands & 
van der Klis 1997]; and SAX J1808.4-3658: i^spin ^ 401 Hz, Az/qpq = 195 ±6 Hz [Chakrabarty et al. 2003; 
Wijnands et al. 2003]). In all other known cases, the largest value of Ai/qpq is approximately consistent 
with either Vspm or i/spin/2. 

3. A mechanism that produces a single sideband is indicated. Most mechanisms that modulate the X-ray 
brightness at two frequencies (such as amplitude modulation) would generate at least two strong sidebands. 
Although weak single and double sidebands have been detected close to the frequency of the lower kilohertz 
QPO (Jonker, Mendez, & van der Klis 2000), at most two strong kilohertz QPOs are observed in a given 
system (van der Klis 2000; Mendez & van der Klis 2000). This suggests that the frequency of one QPO is the 
primary frequency while the other is generated by a single-sideband mechanism. Beat-frequency mechanisms 
naturally produce a single sideband. Because one QPO frequency is almost certainly an orbital frequency, 
the most natural mechanism would be one in which the second frequency is generated by a beat with the 
star's spin frequency or with another orbital frequency. 

4. A successful model should explain why stars with high spin frequencies (as inferred from their 
burst oscillations) tend to have Az/qpq w i'spin/2, whereas stars with low spin frequencies tend to have 

AfQpQ Ri I/gpin- 

All existing beat- frequency models predict that if the bmdamcntal frequencies arc; (1) an orbital fre- 
quency j/orb and (2) the frequency produced by a pattern with n-fold symmetry rotating with the star, then 
the principal beat frequency will be n(forb — t'spin) (Lamb et al. 1985). The sonic-point beat-frequency model 
does predict the appearance of weaker QPOs at a variety of other frequencies (MLP98). One of these could 
by chance differ from the fundamental orbital frequency by about i/spin/2, but this would be an ad hoc ex- 
planation and would not explain a strong QPO at this frequency. Consequently, when evidence arose that in 
many cases fQpQi ~ fQPQ2 — zA)urst/2, it was suggested that perhaps in these stars Z4,urst ~ 2i/spin (MLP98; 
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see also Strohmaycr ct al. 1998b, Muno et al. 2001). The compelling evidence that i^qpoi ~ ^'qp02 — 
in SAX J1808.4— 3658 requires a model in which this frequency arises naturally. At the same time, beat- 
frequency models have many advantages, as discussed above. We propose a beat-frequency model in which 
t'QPOi is either approximately i'qpo2 - I'spin or approximately i'qpo2 - i'spin/2. 

3. The Sonic-Point and Spin-Resonance Model 

We first show that the star's magnetic and radiation fields will excite vertical motion in the disk at the 
"spin-resonance" radius where Vorh — t'spin equals the vertical epicyclic frequency v^, and that the orbital 
frequency at this radius is approximately i'spin/2. We propose that the upper kilohertz QPO is generated 
by clumps orbiting at the sonic point via the same mechanism as in the original sonic-point beat-frequency 
model (MLP98). Indeed, there is new evidence supporting this identification. We then show that interaction 
of the radiation pattern generated by clumps orbiting at the sonic-point radius Vgp with the vertical motion 
of gas at the spin-resonance radius can produce a second QPO at either forb(^sp) — i^spin/2 or Uorhirsp) — I'spin- 
We propose that this second QPO is the lower kilohertz QPO. 

3.1. Generation of motions with frequency i'spin/2 

The magnetic and radiation fields of a neutron star are not perfectly symmetric about its rotation axis. 
Consequently, as the star spins and the gas in the disk orbits, each element of gas experiences periodically 
varying radial and vertical forces. If the frequency of this forcing is close to a natural frequency, the response 
of the element of gas may be large. The two most basic natural frequencies of an element of gas — aside from 
its orbital frequency — are its radial and vertical epicyclic frequencies. Vertical epicyclic motion appears more 
likely to generate a QPO, for two reasons. First, oscillation of an element of gas in the vertical direction can 
generate an X-ray signal by, for example, scattering radiation produced elsewhere in the system into or out 
of the observer's line of sight. In contrast, oscillation of an element of gas in the radial direction does not 
offer such possibilities. Second, vertical motion of gas in an accretion disk may be weakly damped (Markovic 
& Lamb 1998) whereas radial motion is strongly damped (Markovic & Lamb 2000). 

Perturbation of the disk flow by a magnetic or radiation field that rotates with the star will preferentially 
excite vertical motion in the disk at the spin-resonance radius rgr defined implicitly by the resonance condition 

^spin ^orb (^sr) ~ ^tl^iXsr) i (I) 

where i'^(r) is the vertical epicyclic frequency at radius r. Here i/gpin — i'orb(''sr) is the forcing frequency 
experienced by an element of gas orbiting at Vsr- Numerical simulations (see Fig. 1) show that the vertical 
displacement of gas is much greater at the resonant radius than at any other radius. In a Newtonian 1/r 
gravitational potential, v^{r) = Vothii")- In general relativity, i'^(r) is not exactly equal to t'orb(?')) but the 
difference is < 2 Hz at the radii of interest (where forb < 300 Hz). Consequently, at the resonance radius 
where vertical motion is preferentially excited, Vorh ~ Vxp ~ i'spin/2, i.e., at this radius the orbital and vertical 
frequencies are both approximately half the star's spin frequency. 

Each individual element of gas at the spin-resonance radius orbits the star with frequency i'orb(^sr) ~ 
t'spin/2. However, the pattern created by the elements of gas that are above the plane moves around the 
star with frequency ^pattern = fspm- Stated differently, the azimuth where the gas is highest above the disk 
advances at the rate I'spin, even though no gas moves at this rate. The reason for this is that the response of 
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Fig. 1. — Simulated response of gas in the accretion disk to a periodic perturbing force with frequency 
'^spin = 400 Hz. The vertical oscillation of the gas at the resonance radius has frequency fspin/2 and a much 
larger amplitude than at any other radius in the disk. In this simulation the damping time was assumed to 
be 1,000 s. 

the gas is almost in phase with the periodic driving force. For example, if the force on a particular element of 
gas is maximal at the time-dependent azimuthal position 0max(i) = (f'a + t'spin^ and when the gas is furthest 
above the disk plane, then as the forcing field rotates, the gas in the disk at Vgr will be at the peak of its 
vertical excursion at the azimuth (f) = (^max(i)- 

So far we have considered the effect of forces acting only on one (the 'upper" ) surface of the disk. The 
second harmonic of Jospin is weak in the persistent X-ray emission of all known millisecond accretion-powered 
pulsars and, where both have been determined, i^burst is close to I'spin rather than 2i^spin, indicating that we 
observe only one (the "upper") magnetic pole of the star and its X-ray beam as the star rotates. The star's 
opposite pole and its X-ray beam interact with clumps when they are below the lower surface of the disk 
with a phase that will further amplify their vertical motions. 

3.2. Generation of the upper kilohertz QPO 

We propose that the upper kilohertz QPO is related to the orbital motion of clumps at a special radius 
near the neutron star. As shown in MLP98, the radius Vgp of the sonic point in the disk is a strong candidate 
for this special radius. The sonic transition at Vgp is a strong-field general relativistic effect and is usually 
produced by the force on the gas in the disk exerted by radiation from the star. At r^p, the inward radial 
velocity in the accretion disk increases sharply. The sharpness of this change in the radial velocity matches 
the observed narrow width of the upper kilohertz QPO in power density spectra. In addition, radiation 
and magnetic forces are likely to enhance clumping at the sonic radius over clumping further out in the 
disk. As in the original sonic-point model, we assume the frequency i'qpo2 of the upper kilohertz QPO is 
generated by mapping of accretion streams from clumps orbiting at rgp onto the surface of the star. This 
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Fig. 2. Power spectra of the X-ray flux modulation produced by simulations of a disk with a small mimber 
of clumps near the spin-resonance radius (left-hand panel) and a large number of clumps (right-hand panel). 
The star's spin frequency is 400 Hz while the orbital frequency at the sonic point is 1,000 Hz. These power 
spectra demonstrate that if the flow near the spin-resonance radius is clumpy, the effect of individual clumps 
dominates and the dominant frequency is i'orb(''sp) — '^spin/2. If instead the flow is relatively smooth, the 
effect of the clump pattern dominates and the dominant frequency is t'orb(''sp) ~ ^'spin- This simulation did 
not include any signal with the orbital frequency of the gas at the sonic radius. 



mechanism for generating the upper kilohertz QPO is supported by the observed short-term anticorrclation 
of the frequency of the upper kilohertz QPO and the X-ray flux in 4U 1608—52 (Yu, van der Klis, & Jonker 
2001; Yu & van der Klis 2002). 

As discussed in Lamb & Miller (2001), the inward drift of the clumps at r^p causes the footpoints of 
the streams on the stellar surface to revolve around the star with an angular velocity slightly different from 
I'oThifsp)- The radiation pattern they produce rotates with frequency i^Qpo2 = i^ovhifsp) — Si^md^ where Sfrad 
is typically positive and increases in magnitude from nearly zero to a few tens of Hertz as Vgp decreases. 



3.3. Generation of the lower kilohertz QPO 

Scattering of radiation from the neutron star by the clumps orbiting in the disk close to the sonic radius 
will create a radiation pattern on the surface of the disk that rotates with frequency i'orb('^sp)- This pattern 
will preferentially illuminate or shadow the gas orbiting at the spin-resonance radius, when it is above the 
disk plane. The rotating radiation pattern will be scattered much more strongly by clumps orbiting at the 
spin-resonance radius than by clumps elsewhere because clumps at this radius are making large vertical 
excursions, unlike clumps elsewhere in the disk. 

Generation of a lower kilohertz QPO with frequency j^qpoi ~ forb(?'.sp) — i^spin/2.— Suppose first that 
the gas in th(^ disk near the spin-resonance radius is highly clumped. When illuminated, each clump orbiting 
at rsr scatters radiation in all directions. In effect, each clump redirects the radiation propagating outward 
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Table 1: 


Fits to Ai/Qpo vs. 




data 


Source 


M(Mq) 


Az.o(Hz) 


coeff 


xVdof 


Sco X-1 


1.4 


365 


0.0162 


96.3/46 


4U 1608-52 


1.96 


310 


0.0006 


22.1/10 


4U 1728-34 


1.67 


364 


0.0017 


7.3/6 


4U 1820-30 


2.0 


279 


0.00001 


7.2/18 



from the sonic radius in the modest soHd angle that it subtends (as seen from the sonic radius) into all 
directions. From the point of view of a distant observer, each individual clump looks like a light bulb that is 
blinking on and off with a frequency equal to i'orb(''sp) — t'orb(^sr) ~ t'orb(^sp) — t'spin/2. If there are only a 
modest number of clumps at rgr, the effect is somewhat like what one would see if a movie marquee had only 
a few light bulbs, so the blinking of the individual bulbs dominates the time variation. Because the radiation 
is scattered in all directions, an observer does not have to be close to the disk plane to see the X-ray flux 
modulation. The left-hand panel of Figure 2 shows the power spectrum of the flux variation generated in 
a simulation in which five randomly-positioned clumps scatter the radiation pattern coming from the sonic 
radius. The peak at forb(''sp) — t'spin/2 is clearly dominant. 

Generation of a lower kilohertz QPO with frequency uqpoi ~ t'orb(?'sp) — I'spin- — Suppose now that the 
gas in the disk near the spin-resonance radius is less highly clumped. There may be a larger number of 
smaller clumps or the flow may even be relatively smooth. As before, each element of gas is oscillating 
vertically with frequency i'spin/2. Together they form a pattern of raised fluid elements that rotates around 
the star with frequency z^gpin- Because a large number of fluid elements are scattering radiation to the 
observer at any given moment, their individual contributions blend together, so the dominant time variation 
has frequency i^orh{r sp) — Jospin- The effect is somewhat like what one sees when viewing a movie marquee 
with hundreds of light bulbs blinking in phase so that the pattern formed by the lighted bulbs moves around 
the marquee. In this case the brightness variation produced by the pattern of lighted bulbs dominates the 
brightness variation produced by the individual bulbs. The right-hand panel of Figure 2 shows the power 
spectrum of the flux variation generated in a simulation in which 500 randomly-positioned clumps scatter 
the radiation pattern coming from the sonic radius. The peak at Voihi^sp) — ^'spin is clearly dominant. 

Variation of the kilohertz QPO frequency difference. — In the model proposed here, the frequency of the 
upper kilohertz QPO is close to but less than Voih{r sp) and the difference varies by as much as several tens 
of Hertz. However, the frequency of the lower kilohertz QPO is very close to either Voihif sp) — I'spin or 
t^oih{r sp) — Jospin /2, depending on whether the flow at rsp is smooth or clumpy. Consequently, the frequency 
separation between the two QPOs is not exactly I'spin (or i'spin/2). Table 1 shows the results of fits of this 
model to the observed changing frequency separation in several sources (see van der Klis ct al. 1997; Mendez 
et al. 1998a; Mendez et al. 1998b; Ford et al. 1998; Psaltis et al. 1998; Mendez & van der Khs 1999; Jonker, 
Mendez, & van der Klis 2002), adopting the model of Lamb & Miller (2001; note that in that work the 
frequency of the lower kilohertz QPO also varies by a few tens of Hertz, whereas here only the frequency of 
the upper kilohertz QPO varies). The fit is reasonable for all these sources except Sco X-1, which has an 
unusually complicated Aj/qpo-i'qpo2 relation. 
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4. Discussion and Conclusions 

The sonic-point and spin-resonance model proposed liere appears able to explain the observed frequencies 
of the two kilohertz QPOs, the variation of their frequencies by hundreds of Hertz, the closeness of their 
frequency separation Ai/qpo to i^spin in some sources and to I'spin/S in others, and the variation of Auqpq 

with i^Qpo2 (accretion rate) in a given source. 

The difference between "fast" and "slow" rotators. — In the model proposed here, i^qpoi ~ i'orb(''sp) ~ 
^'spin/2 if the resonance with the stellar spin occurs where the gas in the disk is highly clumped, whereas 
I'QPOi ~ ^'orb(''sp) — fspin if the resonancc occurs where the gas in the disk is relatively smooth. Magnetic 
forces may cause the gas in the accretion disk to become more clumped as it approaches the neutron star 
(MLP98; Lamb & Miller 2001). Consequently, the parameters that may be most important in determining 
whether the flow at the spin resonance radius Vgr is clumpy or smooth are the star's spin frequency and 
magnetic field. For a given stellar magnetic field, the flow is likely to be more clumpy if the star is spinning 
rapidly and rgr is therefore close to the star. For a given spin rate, the flow is likely to be more clumpy if 
the star's magnetic field is stronger. 

The four sources with Vgpin > 400 Hz and measurable frequency separations have Ai/qpq ~ Vgpin/^ 
whereas the three sources with Ugpin < 400 Hz have Ai/qpq ~ ^'spin- With such a small sample, one cannot 
make any definite statements, but the apparent trend is consistent with the sonic-point and spin-resonance 
model. The model suggests that if kilohertz QPOs are detected in the recently-discovered 185 Hz, 191 Hz, and 
314 Hz accretion-powered X-ray pulsars XTE J0929-314 (Galloway et al. 2002), XTE J1807-294 (Mark- 
wardt. Smith, & Swank 2003a), and XTE J1814-338 (Markwardt & Swank 2003), their frequency sepa- 
rations should be approximately equal to the respective spin frequencies. The 435 Hz spin frequency of 
XTE J1751— 305 (Markwardt et al. 2002) is high enough that Ai/qpo could be either approximately 435 Hz 
or approximately 217 Hz; QPOs at both frequencies might even be detectable. 

Other oscillation frequencies. — We expect the response at the spin- resonance radius defined by I'orb ^ 
^'spin/2 to be largest, but there may be other, weaker motions at other low-order resonances. In principle, 
weak beat-frequency QPOs could be generated by responses at the radii where forb/j^spin is close to ratios 
of small integers, such as 1/3, 2/3, and so on. The strength of such secondary QPOs would depend on the 
smoothness of the flow and the rates at which vertical motions are damped at these radii. If the flow is 
clumpy, so that the primary beat is with the local orbital frequency, then several beat frequencies could 
be observed. If the flow is instead smooth, so that the primary beat is with the pattern frequency, then 
because the pattern frequency is always the spin frequency, no other significant QPOs would be expected. 
In either case, the stronger the damping is, the weaker the other QPOs would be relative to the one at the 
fundamental spin-resonance radius. 

We note that there is no known reason why the mechanism for producing a lower kilohertz QPO proposed 

in the original sonic-point beat-frequency model would not operate. Apparently this mechanism does not 
produce a strong QPO in the fast rotators, but it might produce a weak QPO in these sources. If it operates 
in the slow rotators, it would produce a QPO near i'orb(?'sp) — ^'spin that might appear as a sideband to the 
lower kilohertz QPO. 

Although the model outlined here is qualitatively consistent with the basic properties of the kilohertz 

QPOs, many aspects of the model require further exploration and development. For example, what if any 
sidebands or additional QPOs are to be expected? How large are the vertical excursions that would be 
expected in a model that includes more of the physics of accretion disks? Is the model quantitatively 
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consistent with the observational upper Umits on X-ray osciUations at Jospin? It wih be important to pursue 
these questions, but in the meantime resonant excitation of gas in the disk appears to be worth exploring as 
a mechanism for generating the lower kilohertz QPO. 

We thank Deepto Chakrabarty and Michiel van der Klis for helpful discussions. This work was supported 
in part by NSF grant AST 0098399 and NASA grant NAG5-12030 at Illinois, and by NSF grant AST 0098436 
at Maryland. 
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